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^ ■ Abstract 



In this paper we reexamine the problem of the separation of spin and charge 



. degrees of freedom in two dimensional strongly correlated systems. We es- 

o 

0> , separation. Specifically, we discuss this issue in the context of the Heisenberg 



tablish a set of sufficient conditions for the occurence of spin and charge 



model for spin-1/2 on a square lattice with nearest (Ji) and next-nearest (J2) 
neighbor antiferromagnetic couplings. Our formulation makes explicit the 



O . existence of a local SU(2) gauge symmetry once the spin-1/2 operators are 

CJ 

replaced by bound states of spinous. The mean-field theory for the spinous 
is solved numerically as a function of the ratio J2/J1 for the so-called s-RVB 
Ansatz. A second order phase transition exists into a novel flux state for 
J2/J1 > {'h/Jijcv We identify the range < J2/J1 < {J2/Ji)cr as the 
s-RVB phase. It is characterized by the existence of a finite gap to the ele- 
mentary excitations (spinous) and the breakdown of all the continuous gauge 
symmetries. An effective continuum theory for the spinous and the gauge 
degrees of freedom is constructed just below the onset of the flux phase. We 
argue that this effective theory is consistent with the deconfinement of the 
spinous carrying the fundamental charge of the gauge group. We contrast 
this result with the study of the one dimensional quantum antiferromagnet 
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within the same approach. We show that in the one dimensional model, the 
spinons of the gauge picture are always confined and thus cannot be iden- 
tified with the gapless spin- 1/2 excitations of the quantum antiferromagnet 
Heisenberg model. 

71.27.+a, 71.30.+h, 74.20.Mn, 75.10.Jm 
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I. INTRODUCTION 



The issue of the possible breakdown of the Fermi-hquid picture is of considerable interest, 
both from a practical and conceptual point of view. The anomalous (with respect to a 
conventional Fermi-liquid scheme) properties of the metallic state of high Tc materials [|l|, 
its proximity to an insulating magnetic state upon doping, might invalidate the use of 
Fermi-liquid theory. The microscopic foundation of Fermi-liquid theory itself has been a 
challenging theoretical problem for a long time. Fermi-liquid theory assumes that there is 
a one to one relationship between the low energy excitations of the non-interacting system 
and the low energy excitations of the interacting system which preserves quantum numbers 
like the spin, charge or particle number. In contrast, it has been proposed that for certain 
strongly interacting systems, the low lying excitations are of a radically different nature 
than in the non-interacting limit. For example, the spectrum just above the ground state 
would be made of "quasiparticles" carrying either the spin or the charge quantum numbers 
of the non-interacting excitations but not both simultaneously. This is the scenario of spin 
and charge separation which has been proposed by analogy to the Hubbard model in one 
dimension by some or to the fractional quantum Hall effect by others [Q. 

The purpose of this paper is to reexamine the conditions under which full separation of 
spin and charge occurs in two dimensional electronic systems. It is believed that a necessary 
condition for this phenomenon to take place is that the ground state in the spin sector is 
a spin liquid P,0J^. A minimal requirement on a state describing a spin liquid is that it 
does not support any long range order breaking the continuuous symmetry under global 
spin rotations, for example Neel ordering. However, under these conditions, many different 
types of spin liquids can still be constructed. Of crucial importance to applications to high 
Tc superconductivity is whether the spin liquid is separated from the excited states by a gap 
or not and in the latter case the nature of the density of states. 

The separation of spin and charge is assumed to occur for a certain range of electronic 
densities as a result of the competition between strong electronic correlations in the form 
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of antiferromagnetic (AF) interactions and the kinetic energy of the electrons. We will be 
interested in the Hubbard model with moderate to strong U or its generalization, the t — J 
model 0], on a two dimensional lattice close to half- filling. At half-filling, the t — J model 
and the U/t —>■ oo limit of the Hubbard model |p reduce to the antiferromagnetic Heisenberg 
model for spin-1/2. The analytical treatment for the separation of spin and charge usually 
starts with the factorization of the creation operator cj^, for the band electrons occuring in 
the Hubbard or t — J model into a pair of operators satisfying fermionic and bosonic algebra, 
respectively |P,|lOl- At this point, the assignment for the spin and charge quantum numbers 



is arbitrary. We will choose the fermionic operator sj^ to carry the spin quantum number 
{a =T; i) but no charge and the state created by is a bare spinon. The bosonic operator 
hi must then carry the charge and the state it creates is similarly identified with a bare 
holon. Our choice is a matter of convenience. However, we believe it is well suited for the 
description of a spin liquid with the shortest possible AF correlations, i.e., of the order of 
the lattice constant. The alternative option consisting of fermionic holes and bosonic spins 
has been used to describe a spin liquid close to the onset of a state with long range Neel 
ordering |TI[ . 



The local phase of the bare spinon and holon operators is arbitrary. Accordingly, there 
exists a local gauge U(l) symmetry once the c|q, = sj^/ij representation is chosen. It thus 
follows that the bare spinon (holon) operators can not create physical states. Indeed, these 
operators break the local symmetry and any theory based on a description in terms of spinous 
and holons must necessarily involve gauge fields. By construction, all of these theories are 
strongly coupled and, thus, the gauge fields vary rapidly. The likely scenario is then the 
confinement of all excitations transforming non-trivially under the gauge symmetry. Thus, 
it would appear that holons and spinous are permanently confined in bound states and that 
there may not be a separation of spin and charge in terms of spinous and holons carrying 
gauge quantum numbers. Nevertheless, several mean-field theories for the bare spinous and 
holons have been constructed at half- filling |T2|-|I3] and away from half- filling |[TP|. The 



challenge to any theory supposed to describe spin and charge separation is to construct 



physical states associated with the spin and charge separately. 

In this paper, we will concentrate on the pure spinon sector at zero-temperature by 
working at half-filling in which case the system is a Mott insulator |20|. A special feature of 



the Heisenberg limit is that the local U(l) symmetry is enlarged to SU(2) pl|j2^ . We will 
refer to this symmetry hereafter as the color symmetry to distinguish it from the symmetry 
under global spin rotations. We also do not have to deal with the complexities caused by the 
dynamical interactions between the spinous and the holons ||23|. We will include frustration 
through a next-nearest-neighbor antiferromagnetic interaction in the Heisenberg limit. The 
issue of spin and charge separation then reduces to the existence of spin-1/2 excitations in 
the low energy sector of the theory. 



Wen ||2J] has proposed a mean-field theory for the strongly interacting spinous which 
preserves the discrete symmetries under time-reversal and parity as well as the space group 
symmetries of the lattice. However, for a certain range of values for the ratio of the compet- 
ing nearest ( Ji) and next-nearest (J2) interactions, his mean-field theory breaks completely 
the continuous part of the local color symmetry. The fundamental assumption is that an 



Anderson-Higgs mechanism |^ , driven purely by the quartic spinon interaction, takes place 



and that its occurrence is sufficient to render the mean-field prediction of the existence of 
spinous in the physical spectrum robust to the strong gauge fiuctuations caused by the con- 
strained nature of the system. This BCS-like ground state is highly unconventional. Pairing 
takes place on the direct lattice and the "spontaneously broken" local gauge symmetry is 
not U(l) but SU(2). It is thus argued that through the macroscopic condensation of spinon 
bound states (the typical size of which is the lattice spacing), the physical spectrum develops 
a gap to all excitations among which some carry the spinon quantum number. Frustration is 
the trigger for this mechanism. A related mechanism has also been proposed by Read and 



Sachdev in their analysis of the Sp(N) Ji — J2 ~ J3 antiferromagnetic Heisenberg model p6 

However, it is not at all clear to us that the spinon/holon picture is the appropriate 
one for the description of spin and charge separation. In fact, to our knowledge, it has 
never been demonstrated that the separation of spin and charge known to occur in many 



one dimensional electronic models belonging to the Luttinger liquids universality class p7 
becomes transparent in the spinon/holon representation. A first example is the Hubbard 
model on a linear chain. There, the separate spin-1/2 and charge sectors are created by 
soliton-like operators which, when expressed in terms of the band electron operator c|q,, 
are highly non-local and therefore cannot be described in any simple way by the spinon 
or holons. Another example is the antiferromagnetic quantum spin-1/2 Heisenberg chain 
with nearest-neighbor interactions for which the lowest excitations above the ground state 
are known to transform like the fundamental representation of SU(2) under spin rotations 
||29| . However, these excitations are again created by soliton-like operators which are related 
to the original spin-1/2 operators through a highly non-linear and non-local mapping, the 
Jordan- Wigner transformation |30|-p^. Nevertheless, the spinon/holon picture offers one of 
the few analytical methods to study the issue of separation of spin and charge provided the 
mean-field predictions are not taken at face value. Moreover, aside from this issue, it offers 
a powerful technique to study non-Fermi liquid behaviour, once the local gauge invariance 
is properly accounted for [P^]. 

In this paper, we will argue that the mechanism thought by Wen can but does not 
necessarily insure that the spinous remain in the physical spectrum. To this end, we have 
constructed an effective field theory for the low lying spinon modes in the presence of smooth 
fiuctuations of the gauge modes. This effective field theory describes the low energy physics 
close to a second order phase transition from Wen's s-RVB phase for < J2IJ1 < {-hlJilcr 
into a gapless phase that, we found, is present at the mean-field level. The gapless state for 



J2/ Ji > {J2/Ji)cv is qualitatively analogous to the fiux state of Affleck and Marston p4 
However, in our problem, the gapless phase occurs at the mean-field level as a result of the 
competition between nearest ( Ji) and next-nearest ( J2) neighbors. Hereafter, we refer to the 



gapless state as the fiux state even though it is not the same state as in [Q. Close to and 
in the fiux phase, the low energy spinon modes and large scale lattice gauge fiuctuations are 
equivalent to a relativistic theory in 2+1 space-time of four Dirac fermions minimally coupled 
to non-Abelian gauge fields. The multiplication of the fermion species (fiavors) is the result 
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of a degeneracy of the mean-field spectrum. Each flavor can be thought of as representing 
low lying fermionic excitations around one of the isolated but degenerate minima of the 
Brillouin Zone. The gauge fields are long wavelength fluctuations of the soft modes residing 
on the next-nearest-neighbor links of the lattice, i.e., the SU(2) color degrees of freedom. 
The nearest-neighbor exchange coupling Ji causes the different fermionic flavors to interact 
through effective scalar flelds. These scalar flelds represent long wavelength fluctuations of 
the nearest-neighbor link degrees of freedom. It is important to stress that these bosonic 
degrees of freedom have no existence of their own independently of the spinous. They 
represent collective modes for the spinous. Just below the onset of the flux phase, the scalar 
gauge flelds acquire expectation values which break all continuous gauge symmetries down 
to a discrete subgroup. Integration of the spinous then yields a local effective theory for 
the non-Abelian gauge flelds and the scalar flelds. We give a non-perturhative argument for 
the existence of the color quantum number of the spinon in the spectrum of the bosonized 
effective theory, provided both kinetic energy scales for the gauge flelds and scalar flelds 
(Higgs) are large enough. The calculation of these two scales will be published elsewhere 
P5| . This mechanism is different from the one constructed from the chiral spin liquid since 
it does not rely on the breaking of parity and time reversal, which at the level of the effective 



continuum action implies the presence of a Chern-Simon term Indeed, although, the 

four Dirac fermions acquire a mass it comes with an alternating sign for each flavor, and 
thus, upon fermionic integration, no net Chern-Simon term results. 

In section |I|, we rapidly describe how to map the Heisenberg model into a SU(2) lattice 
gauge theory. Our strategy to study the issue of spin and charge separation is outlined in 
section |T|. In section |^ the numerical solutions of the saddle-point equations for the s-RVB 
Ansatz proposed by Wen are presented. A generalization of the s-RVB Ansatz to higher 
dimensions is also included. A fleld theory for the low lying excitations in the vicinity 
of the flux phase is derived in section |V|. We contrast our results with an application of 
our formalism to the one dimensional Heisenberg chain in section |V|. Our conclusions are 



presented in section [VII 



II. THE EQUIVALENCE OF THE HEISENBERG MODEL FOR SPIN-1/2 TO A 

SU(2) LATTICE GAUGE THEORY 



The band electron operator of the Hubbard model can always be recast in the form 
^,0 



cni)= |o ><T(i)l +(-)! i(T) ><n I (2.1) 

provided the creation operators , and S for the holon, spinon and doublon states, 
respectively, satisfy the appropriate commutation relations. In the limit of strong on site 
repulsion, the term involving the doublon state is neglected [|T0[. We will choose the spinon 
operator to obey fermionic commutation relations. 

At half- filling, the Hubbard model reduces to the Heisenberg model for spin-1/2 ^j: 

H=Y: ■ (2.2) 

<ij> 

where < ij > is an ordered pair of sites on an arbitrary lattice A, Jij are antiferromagnetic 
coupling constants and 



a being the Pauli matrices. The spinon representation, Eq. ( |2.3| ), for the spin-1/2 degrees 
of freedom must be supplemented with any of the three constraints 

1 = 4a ^ = ^ = S,^ S^^ (2.4) 

for all sites i of the lattice. 

From the fully symmetric tensor dap and the fully antisymmetric tensor eap of SU(2), 
the two bilinear forms 

Xlj = -^L ^jP (2.5) 

and 



Vij ^ia ^af} '^jf3^ 



(2.6) 



can be used to describe a singlet pairing of the two spin-1/2 located on site i and j, respec- 
tively |[15|JT^ . Indeed, the identity 



Si ■ Sj 



1 t 



^ Xij Xij + ^ 1 



1 
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(2.7) 



holds in the Hilbert space of one spinon per site. A spin liquid which, by definition, should 
not show any long range magnetic order, implies, in the spinon picture, the exponential 
decay with separation \i — j\ of the vacuum expectation values < rjjj > or < xlj >■ 

The dynamics of these bilinear forms can be obtained from the vacuum persistence 
amplitude 



Z = Jv [So] Jv [s*]V [s ] e+'Z'^*^' 



(2.8) 



where the lattice Lagrangian is 



L' = ^ita ^dt s 



ita 



J, 



+ E / {VijtUijt + XijtXijt) ■ 



(2.9) 



The integration over the Lagrange multipliers 



(2.10) 



enforces the constraint of single occupancy on the spinon Hilbert space in a redundant 
way. But this redundancy allows for the mapping of Eq. ( |2.9| ) into the Lagrangian of a 
SU(2) lattice gauge theory |]2T]], if a Hubbard-Stratanovich transformation with respect to 
the composite fields ri and x is performed first. 
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With a particle-hole transformation of the spinons, our final Lagrangian will then take 
the form ^ 

L = J2^u {idt-A,u) At (2.11) 

i 

<ij> 



Here, the V^'s are 



At=\ . (2.12) 



The Aa's belong to the fundamental representation of the su(2) Lie algebra 

Aoit = ^ Sou ■ <?. (2.13) 



Finally, the W^s are 2x2 matrices of the form 

W,,,= \ J 1, (2.14) 



which satisfy 



rp* I v* 



W,,, = Wl,. (2.15) 



The entries E and X of the W's, are the Hubbard-Stratanovich degrees of freedom associated 
with the spinon bilinears ?7 and Xi respectively. 

The lattice Lagrangian in Eq. ( p.ll| ) is left unchanged by the local gauge transformations 

V'it ^ Ait = UitAt, 

A,,, ^ a;,, = U,, A,u Ui + iidt f/J Ul 

W,^, - Wl^, = U,, W,^, ul (2.16) 

for all f/jj G SU(2). This local symmetry will be called a color symmetry. It is a different 
symmetry from the one generated by global spin rotations. Indeed, under the particle-hole 
transformation Eq. ( p.l2| ), the spin-1/2 operators of Eq. ( |2.3| ) are mapped into 
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-^f = +1 -1)^+1 (m, -1). (2.17) 

The bilinears b and m defined above are left unclianged by tlie local color transformation 
Eq. (|2.16|) and thus the Heisenberg interaction explicitly transforms like a color singlet when 



expressed in terms of the ip^s. Spin-spin correlations can be obtained from the generating 
functional 

Z[J] = Jv [W^] V [W]Jv [as] J V [r]V[^ ] 

xe+'K ^ + S« ^"^-^^ ) (2.18) 

where the source term is written in terms of the bilinears b and m defined by Eq. ( |2.17| ). 

It is important to realize that the gauge degrees of freedom in Eq. (|2.11| ) are not indepen- 
dent from the fermionic degrees of freedom. For example, neither A^, nor the SU(2) factor 
of Wij possess a lattice version of the kinetic energy. Their presence is simply a device to 
project the Fock space which is generated cyclically from the vacuum state |0 >^ defined 
by 

V'm |0 >^ = 0, Vz e A, a = 1,2, 

onto the physical Hilbert space which is the tensorial product over all sites i of the vector 
spaces spanned by the color singlet states |0 and 6||0 (see p7| ). 



III. SPIN LIQUID, DECONFINEMENT AND SPIN AND CHARGE 

SEPARATION. 

In the previous section we have made use of the mapping of the spin-1/2 antiferromag- 
netic Heisenberg model into a SU(2) lattice gauge theory with the hope that the latter 
formulation is better suited for the description of a spin-liquid ground state and the possible 
existence of spin-1/2 excitations. The fact that the problem is, in effect, a gauge theory 
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indicates that in all likelyhood the spinons will necessarily be permanently confined inside 
SU(2) color singlet bound states. This is so since, below four space-time dimensions all 
(compact) gauge theories are generally in a confined phase. If this is the case, there is little 
hope that this mapping may be of great help in producing a mechanism for the separation 
of spin and charge quantum numbers in the physical spectrum of the system. Hence, the 
problem that needs to be solved is to find mechanisms that separate the spin and charge 
either in a manner consistent with the confinement of spinons or by dynamically driving 
the system into a deconfined state. In this section we will show that Wen's s-RVB state 



effectively realizes the second scenario. In section 0, we show that the first scenario is 
realized in the one dimensional version of Wen's state. 

In order to describe a spin liquid, an eigenstate of the Hamiltonian of Eq. (|2.2|) should 
show exponential decay with separation for all spin-spin correlations which break the global 
rotational spin symmetry. For example, to avoid Neel ordering we can imagine partitioning 
the lattice into distinct ordered pairs and choosing for each such pair the singlet state of the 
four dimensional Hilbert space associated with the two spin-1/2 degrees of freedom. The 
expectation value of Si ■ Sj in this state is only non-vanishing when i and j belong to the 
same ordered pair. However, such a state, the so-called valence bond state [^], breaks the 
lattice symmetry. If we require "featurelessness" from a spin liquid, then an easy remedy for 
a valence bond state with nearest-neighbor pairing, say on a square lattice, is to demand a 
uniform expectation value of the valence bond ordering operators for any spin and its four 
nearest-neighbors. This leads us naturally to the RVB overcomplete basis for Hamiltonian 
Eq. ( p.2|) . However, the RVB basis is extremely difficult to work with. 

On the other hand, it appears to be very easy to implement featurelessness and expo- 
nential decay of the spin-spin correlations in the gauge picture, Eq. ( |2.11[ ). Featurelessness 
amounts to selecting only those gauge equivalent classes of configurations C{^.^ w} which are 
closed under the action of the lattice point group, i.e., gauge configurations on which the 
action of a point group transformation is equivalent to a gauge transformation. Exponential 
decay of the spin-spin correlations requires that within the class w} 
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\det Wij\ oc e-l^-^'l (3.1) 

is satisfied, although this condition might not be sufficient (overcompleteness of the RVB 
basis) . 

This suggests the following mean-field approximation for a spin liquid: 

(i) Choose a gauge field configuration consistent with the requirements of a spin liquid. 

(ii) Solve the fermionic sector in the background of the gauge fields given in step (i). 

{Hi) Impose the self-consistency condition that the background gauge fields given in step 
(i) are the expectation value of the appropriate fermionic bilinear forms in the ground state 
found in step (ii). 

The drawback of this approach is that the self-consistency condition only satisfies the 



local constraint of single occupancy for the spinous, Eq. ( p.4|) , on average. In the worst 
case scenario, one could imagine that all the eigenstates of the mean-field Hamiltonian for 
the spinous belong to the unphysical Hilbert space, although on average the constraint of 
single occupancy is satisfied. For example, at half-filling, a state with all but two sites singly 
occupied by spinous is orthogonal to the physical subspace of the fermionic Hilbert space. 
In the language of the ip^s, single occupancy on any given site is the local constraint of 
Gauss's law required by the local gauge invariance, namely the local physical states are the 
color singlets |0 >^ or 6| |0 >^,. It is then easy to convince oneself that the ground state 
and the single particle excitations of the mean-field theory have always a non-vanishing 
projection onto the unphysical Hilbert space. Indeed, the mean-field ground state is a Slater 
determinant for the ^'s which can be obtained by a Fourier transformation of 

n |0 (3.2) 

where Qq denotes the set af all |A|/2 wave vectors entering in the Slater determinant of the 
mean-field ground state. The mean-field ground state and the excited states are thus con- 
volutions in direct space whereas the Fourier transformation of the gauge invariant operator 
bj = 4'liipi2 is a convolution in reciprocal space. 



13 



Since the problem of the quantum numbers of the physical states is non-perturbative by 
nature, it is important to take into account the quantum fluctuations of the gauge fields 
around the mean-field. The problem we end up with is that of fermions which are strongly 
coupled with gauge fields. If the mean-field theory predicts a gap in the fermionic spectrum 
a simpler picture can be drawn. In this case it is possible to integrate out the fermions to 
obtain a local theory. Although the resulting effective theory will be constructed from purely 
bosonic dynamical fields (which can be viewed as local fermionic composites), its spectrum 
can contain states with fermionic quantum numbers. We will show below that this is the 
case for the s-RVB Ansatz where the effective action involves gauge fields locally coupled 
to bosonic matter (scalars) which represent bound states of the spinous. We will also see 
that these bound states have simple transformation properties under the symmetries of the 
system. Of particular importance is that we will be able to find bound states that transform 
non-trivially under the local SU(2) invariance. Hence, at least qualitatively, our problem 
is closely related to that of finding the nature of the ground state of scalar fields coupled 
to SU(2) gauge fields. This is a problem that has been investigated quite extensively in 
the context of Lattice Gauge Theory. In what follows we will review those results, make 
connections with our problem and then draw conclusions. It will be apparent from our 
discussion that the key to the mechanism is to determine the transformation properties of 
the bound states. 

Let us set the notation and define a set of scalar fields $j on the sites i oi a, d + 1 
dimensional space-time hypercubic lattice A and a set of SU(2) gauge fields Vij on the 
ordered links < ij > of the same lattice. The gauge fields transform like group elements 
whereas the matter fields will transform under some representation of SU(2). In practice we 
will consider two ) matter in the fundamental (or spinor) representation of SU(2) 

color and b) matter in the adjoint (or vector) representation of SU(2) color. A simple but 
generic action for this system is 




<ijkl> 
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<'j> 



i i 

p>0, K > 0, A > 0. (3.3) 



The first term of Eq. (|3.3| ), wliicli is a sum over oriented plaquettes < ijkl >, represents 



tlie gauge invariant kinetic energy for the SU(2) gauge fields while the rest of the action 
describes the dynamics of the matter and their coupling to the gauge fields. The scalar 
product of the matter fields $ in a given representation is denoted by (■,■). Finally, we have 
also allowed for a fluctuation of the norm of the matter field which, however, is irrelevant 
to the problem of confinement. 

The properties of the phase diagram in the /? — /t — A space for Eq. ( |3.3| ) which are 
important to us depend only on the dimensionality of space-time and on the transformation 
properties of the matter fields [^,^]. We will distinguish two situations. In one case, the 
matter field $ is in the fundamental representation, in the other it is in the adjoint. Those 
two representations differ in one crucial aspect. The fundamental representation is one to 
one whereas the adjoint is two to one since it fails to distinguish between the two elements 
of the center Csu(2) of the group SU(2) color . 



We first start our exploration of the phase diagram from the point (/?, k, A) = (0, 0, 00). 
By letting A — > 00, we have frozen the norm of the matter field to the value one. When we 
vary (3 from to 00, we obtain a pure, non-Abelian lattice gauge theory. Along this line, 
the Wilson loop Vr^(a) for sources in the fundamental (adjoint) representation 

^V) = < n V^J > (3-4) 

<ii>er 

where F is a closed path on the hypercubic lattice, can be used to detect a phase transition. 
For an appropriate choice of F, 

F = al (3.5) 

the potential 
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V^f(a) (R) ^ - lim hnW f(a) (T, R) (3.6) 

1 — ^oo ± 

measures the energy of two static (infinitely heavy) sources in the fundamental (adjoint) 
representation which are separated by a distance R. 

For small values of P, i.e., strong gauge coupling, the static potential is linear (accordingly 
the Wilson loop obeys an area law) and the static sources are confined. The lower critical 
space-time dimension is 4 for non-Abelian groups. We are interested in space-time of dimen- 
sion 2+1 and therefore do not encounter a phase transition suggesting that for larger values 



of [3 deconfinement (the Wilson loop obeys a perimeter law) of the static sources occurs 

At the point (oo,0,oo), the pure gauge theory freezes into any of the classical gauge 
configurations with all plaquette variables set to unity. The simplest choice is to have all 
Vij set to unity. We then start moving along the line {oo, k, oo). The action for the matter 
field in the fundamental (adjoint) representation and with frozen norm is equivalent to that 
of a classical 0(4) (0(3)) Heisenberg model. For space-time dimension larger than 2 and 
irrespectively of the matter field representation, a second order phase transition will occur 
when K is larger than k^^^^ |^ . 



As of now, the choice for the representation of the matter field had no impact on the 
existence of a phase transition. However, along the line (/?, oo,oo), the physics will differ 
dramatically. At (oo, oo, oo), the matter field in the fundamental (adjoint) representation are 
frozen into a classical configuration, say $0) which minimizes the action and breaks the global 
0(4) (0(3)) symmetry, provided + 1 > 2. By decreasing the value of (3, the gauge fields 
are unfrozen. However, the matter field configuration $o acts like an external symmetry 
breaking magnetic field through the hopping term. Indeed, under a gauge transformation 

($0, V.,<f o) ^ ($0, U,V,,Ul'%). (3.7) 

The adjoint representation differs from the fundamental representation in that if one 
chooses to multiply Vij by the representative of — a" on the left and +0"" on the right, then 
the hopping term remains unchanged in the adjoint while it changes sign in the fundamental. 
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In other words, there will always be a residual local Z2 symmetry along 00, 00) in the 
adjoint representation. Hence, if the matter field breaks all the continuous symmetry, no 
phase transition takes place in the fundamental representation, while a first order phase 
transition takes place at in the adjoint representation |^ . 



Finally, along the line (0, k, 00), the action describes a theory of non-interacting links 
(see [^) which does not undergo a phase transition. 



We can extract valuable informations from Eq. (|3.3|) in the limiting cases where either 
the matter field or the gauge degrees of freedom are passive bystanders to the dynamics. 
In the interior of the phase diagram of Fig. |I], gauge and matter degrees of freedom do not 
decouple from each other. Consequently, the Wilson loop for static charges carrying the 



same representation as the matter field always obeys the perimeter law [46]. Indeed, pair 



creation of matter, which causes the screening of the static color charges |^^, are induced 
by quantum fiuctuations. 

It also follows from this argument that for the matter in the fundamental representation, 
all the Wilson loops, irrespectively of the representation of the static sources, will obey the 



perimeter law in the interior of the phase diagram Wilson loops cannot then be used 



anymore as order parameters. In fact, it was proven that there exists a strip of 

analyticity joining the small (3 and small k region {Confinement region) to the large (3 and 
large k region [Higgs region) of the phase diagram, in which all the gauge invariant Green 
functions are analytic. In short, there is only one phase, the so-called Confinement-Higgs 
phase, in the phase diagram (/?, k, 00) when the matter is in the fundamental representation. 
A line of first order phase transitions emerges from the point (oo,k|!,oo). But as in the 
Temperature-Pressure phase diagram for hquid and gas, it ends up at a critical point in the 
interior of the (3 — k plane (see Fig. |l](a)). 

This analysis implies that if the semiclassical expansion around a mean-field Ansatz 
yields an effective bosonic action with the matter in the fundamental representation, there 
will not be any states in the spectrum which carry the spinon quantum numbers. Conversely, 
if the matter is in the adjoint representation and breaks all the continuous symmetries, then 
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a phase in which there are states that carry the fundamental color and spin-1/2 quantum 
numbers is possible. Indeed, the Wilson loop in the fundamental representation does not 
need to follow the perimeter law everywhere in the interior of the phase diagram. The 
dynamical matter field induced by quantum fluctuations carry the adjoint color charge and 
therefore cannot completely screen the fundamental charge of the static sources. Gauge 
invariant Green function are separately analytic in the confinement and Higgs region. On 
the upper boundary of the Higgs region, the Wilson loop obeys the perimeter law of the 
deconfining regime of a pure gauge Z2 theory. By continuity, the Wilson loop must obey the 
perimeter law in the Higgs region. The same argument applies to the confinement region 
where now the Wilson loop satisfies the area law. There must be a line of phase transition 
separating the two regimes and connecting the classical 0(3) phase transition to the pure 
gauge Z2 phase transition [^] (see Fig. |l|(b)). 

When we allow the norm of the matter field to fluctuate, i.e., A < 00, changes occur in 
the (3 — K, phase diagram. However, the deconfinement and Higgs regimes always remain 
analytically connected for matter fields in the fundamental representation. Conversely, the 
phase transition survives in the adjoint representation pU|. 



This is the important insight that the gauge picture provides us with. The problem we 
are now faced with is how to justify the use of action Eq. ( p.3|) where the matter is in the 
adjoint representations on the basis of Eq. (|2.11|) . This we start now. 



Let us reexamine Eq. (p.ll|) by parametrizing the W^s according to 

W,^, = +i ( v'idet W,^,\ ) V^,^,, V;^, e SU(2). (3.8) 

The local gauge symmetry only affects the SU(2) factor of the link degrees of freedom. As 
it stands, Eq. ( p.ll| ) does not contain a kinetic term for the Ag's and V^s. This corresponds 



to the infinitely strong gauge coupling limit, /5 = 0. This is to be expected from gauge fields 
implementing a local constraint. 

So far, we have not made any approximations. However, in some mean-field states (in 
particular in Wen's s-RVB) there are interesting operators which acquire expectation values. 



For instance, a term which is absent from Eq. but which would be allowed by the 

symmetry, is 

-M.,t = ^ tr ( P,, ly,^, 4 W}^, ) , (3.9) 

for any ordered pair < i j > of lattice sites. The main point of Wen's Ansatz is that the 
operator P^^, the path ordered product 

n 

= n W^.-rk,, lo^l, In^l, (3.10) 

1=1 

(along the ordered pairs < io h >,■■■,< k-i k >,■■■,< in~i in >) does acquire an 
expectation value in the s-RVB ground state. Be aware that the TVI's do not depend on the 
determinant of the VT's. If we take i and j to be nearest-neighbor with e being the lattice 
constant, say j = i + e expand Pjt and 

ly,^, = +i (v^ldet W l) e+^^^«s e su(2), (3.11) 

in powers of the lattice constant, then 

i {M,^, + M*^,) = 1 - ^ tr {{D, P) (D, P)t) 

+ 0(e3), (3.12) 

where is the covariant derivative in the adjoint representation, i.e., 

= d, + i[A,, ■]. (3.13) 

The P's can therefore be interpreted as matter fields in the adjoint representation of SU(2) 
color. The fluctuations around the s-RVB ground state should contain terms of the form 
of Eq. (|3.12|) . Thus, we argue that the fluctuations around an s-RVB state should behave 
qualitatively like a system of matter fields in the adjoint (vector) representation of SU(2) 
(which break all continuous symmetries) coupled to SU(2) gauge fields. We are then un- 
der the general conditions of the theorem relating the confinement and the Anderson-Higgs 
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mechanisms of Fradkin and Shenker Notice, however, that the parameters of the effec- 
tive Lagrangian need to be determined from a microscopic calculation. These parameters 
will, in turn, determine in which phase the symmetry is realized. 

In summary, we have argued how deconfinement could occur for the fundamental charges 
(the spinous) in the phase diagram of a generic lattice gauge theory coupled to matter fields 
in the adjoint representation of the gauge group. Although Eq. (|2.11|) does not resemble a 
canonical SU(2) gauge-matter lattice theory, say for example Eq. ( ^73|) , we have identified 
candidates for a matter field namely the P's. Our goal in the next section will be to 
construct a mean-field theory for a spin liquid such that P's are generated which break the 
gauge symmetry down to Z2 



IV. MEAN-FIELD THEORY. 



A. The s-RVB Ansatz. 



From now on, we will restrict ourself to a square lattice (c?=2) with nearest ( Ji) and next- 
nearest (J2) neighbor antiferromagnetic interactions. The s-RVB Ansatz assumes that the 
true ground state of the system described by Eq. ( p.ll|) is characterized by the configuration 



ijt 



(4.1) 



-Xa^ a j = i + X, 

-Xa^ ifj = i + y, 

-Re E + lm E a"^ if j = i + x^, 
_ —Re E — Im E a'^ if j = i + x_ , 
of the gauge fields [^. The lattice spacing is taken to be one and x = (1,0), y = (0, 1), 
x+ = (+1, +1) and x^ = (+1, —1). The mean-field spectrum ±|^| for the s-RVB Ansatz is 
then given by 



al 



= —-^ + J2 Re E cosfc^cos/Cy, 



~ + J2 Im E sin kx sin ky, 



(4.2) 
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il = + JiX - {cosK + cos ky). 

Because of the particle-hole symmetry, the mean-field ground state is the state obtained 
by filling up all the one-particle eigenstates in the Brillouin Zone Q with negative energy 
eigenvalues: 

I ^s-RVB > = n A ^2 I >^ (4.3) 

k, \U<o 

k 

where | >^ is the state annihilated by all the ip^s. 

As shown by Wen, all the discrete symmetries of the lattice as well as parity and time 
reversal are preserved by the s-RVB Ansatz when a? = a? = 0. But, most importantly to 
us is that the Ansatz is designed to break all the continuous color symmetries provided the 
complex mean-field parameter E is neither real nor purely imaginary and X is non- vanishing. 
Indeed, the counterclockwise (ccw) and clockwise (cw) products 

pccw _ Tir \ir Tir 

^ it — ^^i{i+x)t (i+x)(i+x+y)t (i+x+y)if> 

p™ = ^(i+x){i+x~-y)t ^{i-x+y)ity (4-4) 

reduce to 

pccw _ j^j^g E - Im E a^) , 

pew ^ i^j^g E a^ + lmE a^) , (4.5) 

in the mean- field. Hence, P^'^'" and P™ do not commute if and only if X, Re E' and Im E 
are all non- vanishing. 

The s-RVB is the phase in which X, lie E and Im E are non- vanishing. In this phase, 
|.^^| is also non-vanishing over the entire Brillouin Zone. The spinous have acquired a gap. 
The spin-spin correlation function can be calculated, within the mean-field approximation. 



following the methods of Marston and Affleck [Q. The spin-spin correlation function has 
an exponential decay (in imaginary time) as a function of separation since the s-RVB state 
has a gap. Also, as expected, in the BZA (J2/J1 — > 0) regime the structure factor has a 
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peak at (tt, tt) and, in the flux regime (J2/J1 > {J2/ J\)ct) it has peaks at (7r/2,7r/2), (0, tt) 
and (tt, 0). 



B. The saddle-point equations and mean-field ground states. 

The saddle-point equations in the inflnite volume limit are 

= F,(M), i = l,---,6, (4.6) 
where {n} is the set made of the mean- field parameters 

Ki = ai, K2 = a?, K2, = al, 

K,4 = Re E, K,5 = Im E, hq = X, (4.7) 

and the functions Fj are defined implicitly by integrals over a square with the area 47r^ (the 
Brillouin Zone f2): 

Fv- 
F2-- 



Fs- -J-, !^, F, = — l^g, - K,. (4.8) 



Here, 



' 1 

ATT^Jn 1^1' 


47r2 X 


6 


— fi;4 


1 f 6 


~ 47r2 X 


— 92 


- 


1 f 6 


~ 47r2 L 




- 


9i(k) = 


COS cos ky ^ 






92{k) = 


sin kx sin /cj^. 






93{k) = 


- (cos/Cj; -|- cos ky), 





(4.9) 

each transform irreducibly under the point group of the square lattice. Hence, the integrals 

— * 

only need to be performed over the reduced Brillouin zone Q' — {k e ^ \ < \ky\ < n}. 
It is also useful to consider the Jacobian 

dF 

^'''^d^ ^,i = l,---,6. (4.10) 
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There are two limiting cases for which the saddle-point equations can be solved analyti- 
cally. The BZA limit is the limit when J2 = ||12|. The mean- field parameters 



Ki = K2 = K,3 = K4 = = 0, Kq 



(4.11) 



are (up to a sign) the unique solutions of the saddle-point equations. The single particle 
excitation spectrum is that of a tight-binding fermionic gas at half-filling. The edges of the 
square with vertices (vr, 0), (0,7r), (— vr, 0), (0, — vr) at which the spectrum is gapless is the 
Fermi surface. The Jacobian is formally given by the matrix 

/FiiOOO 0\ 








V 



F2,2 










-^5,5 










F4,4 



Fq^qJ 
It turns out that the non-vanishing entries are divergent integrals. 

The limiting case Ji = can also be solved analytically. The mean-field parameters 



(4.12) 



BZA 



K4 



1^2 



K3 = Kg 
1 1 



0, 



47r^ 



7; I ^9! + 9I ~ 0.339, 
2 Jn * 



(4.13) 



are (up to signs) the unique solutions of the saddle-point equations. The single particle 
excitation spectrum has four discrete gapless points on the Brillouin Zone. They are located 

at 



ko3 



02 



71 

"2 ' )' ^04 



TT 

TT 



(4.14) 



In the neighborhood of each of these nodes, the excitation spectrum takes the same form as 
that of a relativistic massless particle. It is shown in appendix ^ that this gapless state is 
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-^2,2 

















-^3,3 

















-^4,4 


-^4,5 











-^5,4 


-^5,5 


I 















(4.15) 



characterized by an average mean-field flux of vr through a cell bounded by the next-nearest- 
neighbor links and hence it will be called the flux state [^. In this limit, the Jacobian is 
well defined and given by 

\ 





^6,6 /flux 

For arbitrary values of the ratio J2/ Ji-, the saddle-point equations can only be solved 
numerically. The value of the solutions of the saddle-point equations as a function of J2I J\ 
are plotted in Fig. ^ One observes that ^2 = ^3 = 0. This results from the lattice symmetry. 
It also appears that there exists a critical value for the ratio J^j J\ above which the flux mean- 
field solution is always preferred. The numerics predict this critical value to be between 1.3 
and 1.4. 

This critical value can be found analytically by requiring the Jacobian determinant to 
satisfy 

It is understood here that the Jacobian is calculated for the mean-field parameters in the 
flux limit but for arbitrary values of the ratio J2I J\- This condition is thus simply 

J-i 1 1 



9l 



1.342 



(4.17) 



if K4 ^ 0.339. 

The numerics close to the BZA regime are too unreliable to predict the existence of a 
critical value of the ratio J2/ Ji below which the BZA state prevails. The analytical argument 
used for the flux phase fails as well in this limit since the Jacobian is not of maximal rank 
and moreover is ill defined. Indeed, singularities associated with the Fermi surface renders 
all non-vanishing entries of the Jacobian infinite. 
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The numerical results clearly indicate that, between J2/J1 = 0.6 and J2I J\ = 0.7, the 
topography of the single particle energy crosses over from a regime with the absolute minima 
located on the diagonal of Q' to a regime with the absolute minima located on the edges of 
the reduced Brillouin Zone (see Fig. |^). There should therefore exists a crossover value for 
the ratio J2/J1 for which there is an accidental degeneracy of the minima. The Lagrange 
multiplier \al\ appears to keep track of this crossover, the point at which it peaks. Our 
mean-field results are summarized in Fig. 

C. Higher dimensional generalization. 

We will try in this section to generalize the s-RVB Ansatz to any spatial dimension 
d > 2. Let A = Z'^ be a d dimensional hypercubic lattice on which Eq. ( p. Ill) is defined. It 
is helpful to introduce the 2d unit vectors 

x^^ = (0,---, 0, a^, O,---, 0), 

/i = l,---,rf, a^ = ±l (4.18) 
and the 2d{d — 1) vectors of norm \/2 

xj^;^ = (0,---, 0, a^, O,---, 0, a,, O,---, 0), 

H,u = l,---,d, cr^^ = af,-a^, n < (4.19) 

which connect any given site with all its nearest-neighbors and next-nearest-neighbors, re- 
spectively. The straightforward generalization of Eq. ( f4.1| ) is 



-Xa^ ii j = i + xit , 

-Re E a^ + lmE ii j = i + x^\, (4.20) 

^ -Re E a^-lmE ii j = i + x^- 
The breaking of all the continuous color symmetries is automatically achieved once this 

is so for d = 2. However, we need to check that the Ansatz is, up to gauge transformations, 

invariant under all point group transformations of A as well as under time reversal. 

25 



Time reversal is equivalent to complex conjugation. In turn, complex conjugation is 
equivalent to Im E — Imii^. But, this change can be undone by the gauge transformation 

U, = sgn(z) a\ (4.21) 

provided Ao = ^ al cr^. Here, sgn(i) is the lattice function defined in Eq. ( [4.24|) . 

We now turn to the point group transformations of A. They are generated cyclically 
from the d\ permutations of the coordinates of i G A: 

(^1 ! ■ ■ ■ ) ^d) ^ (^ni ) ■ ■ ■ ) ^n^j); 

ni,---,nrf = l,---,ci, (4.22) 

and from the d transformations 

(ii , ■■■ , id) ^ (-«i , ■ ■ ■ , id), ■ ■ ■ 

■ ■ ■ , ■ ■ ■ , id) ^ {ii , ■ ■ ■ , -id)- (4.23) 

When d = 2, the d\ = 2 permutations are the identity and the reflection about the 



diagonal x = y. They leave Eq. ([4.20|) unchanged. The remaining point group transfor- 
mations are generated by reflections about the x or y axis. For d = 2, they are equivalent 
to complex conjugation and Wen's assertion that Eq. ([4.2CI| ) describes a spin liquid follows. 
But for higher dimensions, the equivalence between complex conjugation and any of the d 
reflections about a d — 1 dimensional hyperplane, Eq. ( ^.23| ), is not true anymore. 



A possible cure to this problem is to find a partitioning of all the next-nearest-neighbor 
links into two disjoint families which are invariant under all point group transformations. 
Notice first that the lattice function 

sgn : A ^ Z2, i sgn(i) = (-1)^1+-+*^ (4.24) 

is unchanged by any permutation of the coordinates of i. Furthermore, sgn(i) is also left 
unchanged when any one of the coordinates of i has its sign changed. Consequently, the 
partitioning of A into the even (odd) sites A^j^^^-, is not affected by the point group transfor- 
mations. 
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We immediately see that the Ansatz 



A 



oit 



A, 



o' 



(4.25) 



if J 



i + {ft 



ijt 



-Re Ea^ + lm E a'^ if j = z + x'^^, i G Ae, 
^ -Re Ea^ -ImE if j = i + x'^^, i e Ao, 
respects all the point group symmetries. On the other hand, translational invariance and 
time reversal appear to be broken. However, they affect the Ansatz Eq. ( [4.25| ) in the same 
way as the gauge transformation Eq. ( [4. 2 ID does. We conclude that this Ansatz describes a 
spin liquid which breaks all the continuous symmetries of color SU(2). 

In two spatial dimensions, this new Ansatz reduces to the BZA Ansatz in both the 
J2/J1 and J1/J2 limits. Moreover, it is not possible to construct a mean- field 
plaquette term threaded by a flux of vr. We therefore expect that a gap to the mean- field 
spinon excitations is generated for all < J1/J2 < 00. 



V. GAUGE FIELDS FLUCTUATIONS CLOSE TO AND IN THE FLUX PHASE. 

It is fortunate that the mean-field theory yields the flux phase since one can relatively 
easily identify a continuum effective theory for the long wavelength fluctuations of the soft 
modes. In contrast, it is much harder to do the same close to the BZA limit. This is so 
because in the flux phase the locus of points in the Brillouin Zone for which the mean- 
field single particle energy reaches a minima are isolated. From the location of the isolated 
minima, one can then extract a reciprocal vector around which the Fourier transform of 
the spinon and link degrees of freedom can be expanded, thus yielding the intermediate 
continuum limit. However, it is not feasible to carry out such an approach at the BZA point 
since there there are infinitely many degenerate minima along a square, the spinon Fermi 
surface, with vertices (±7r, 0) and (0, ivr) in the Brillouin Zone. Moreover, the difficulty is 
compounded by the existence of saddle points in the energy surface at the vertices of the 
spinon Fermi surface causing Van-Hove singularities in the mean-field density of states and 
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by the nesting property of the Fermi surface. One possible drastic approximation is to replace 
the anisotropic Fermi surface altogether by a circle and then expand the nearest-neighbor 
links degrees of freedom in powers of the gauge fields up to second order. The continuum 
effective theory is that of gauge fields coupled to non-relativistic fermions. But, by doing so, 
one neglects magnetic instabilities induced by the nesting property of the Fermi surface at 
half-filling. Away from half-filling, such an approximation has been used by Lee and Nagaosa 
33| in their study of the t — J model. It is then known that the fluctuations of the gauge 



fields, when coupled to the "Fermi surface", lead to severe infrared divergences which 
possibly signal an instability of the mean-field theory (and thus of the "Fermi surface" ) to a 
"strange metal" [^. Given these difficulties, we will concentrate on the flux phase and will 
not discuss the BZA regime any further. In any event, the mean-field theory predicts that 
the BZA and flux regimes are smoothly connected. Thus, we expect that a detailed study of 
the flux regime can throw light on the behavior near the BZA limit as well. Notice, however, 
that the low energy fermionic states have different physical properties in both regimes given 
the level crossing of the excited states shown in the previous section. Nevertheless, since the 
mean-field ground state evolves smoothly (i.e., without level crossings) we will still assume 
that the physics is qualitatively similar at both ends. 

In the spirit of section p| , we need to identify the bosonic degrees of freedom of the 
continuum limit. To this end we will perform a semi classical expansion of the link degrees 
of freedom around the mean-field Ansatz. In the language of critical phenomena, we will 
use the fact that the mean-field theory of the previous section predicts a second order 
phase transition at {J2/ Ji)cv We will then think of the continuum limit as the tuning of 
the coupling J2/J1 toward the point {J2/Ji)cv where a correlation length (the inverse of 
the spinon mass) becomes infinite. In this scheme, the intermediate effective field theory 
we obtain is the (infrared unstable) fixed point action of Eq. ( |2.11| ) in the sense of the 



renormalization group (see Fig. Notice that the semiclassical expansion is not fully 
justified since the theory, as it stands, does not have a small parameter that will control 
the fluctuations. Formally, we can regard our semiclassical expansion as a expansion 
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in which has been set to unity. While this approach seems tempting, we should be aware 
that at small values of the physics may be different from the semiclassical (mean-field) 
picture. 

The intermediate effective theory so obtained will contain gauge fields and matter fields, 
which both transform like the adjoint (vector) representation of SU(2) color, and are coupled 
with (several species of) Dirac fermions. However, the bosons are still not truly dynamical 
since they do not possess a kinetic energy of their own. On the other hand, all the sym- 
metry requirements that we are after are fulfilled, namely the bosonic fields break all the 
continuous gauge symmetries down to a discrete subgroup Z2. In the next subsections we 
derive the effective continuum theory in the vicinity of the flux phase. It turns out that 
Wen's choice of gauge is not very natural in the flux regime. In appendix ^ we give a more 
convenient construction in a different ("flux") gauge which makes simpler the description 
of the fluctuations around this state. The reader is referred to appendix ^ for the notation 
that we will use below. 

A. Continuum limit in the flux phase when J\ = 0. 

We will assume that J\ has been turned off. The fermions on the even sublattice are then 
decoupled from those on the odd sublattice. We shall show how to take the continuum limit 
so as to obtain relativistic Dirac spinors in 2+1 dimensions which couple via the minimal 
coupling to two independent set of SU(2) gauge fields. 

We will choose to parametrize small fluctuations of the Q's (see Eq. ( |A12| )) around the 
mean-field Ansatz by 
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(5.1) 



Here, e'(e) is related to the lattice constant e by 



e' = V^e, ( e = 2 e' ). 



(5.2) 



Only SU(2) color fluctuations of the link degrees of freedom are considered since the deter- 
minant l^'l is kept unchanged. 

In the Heisenberg representation, the equations of motion satisfied by the fermions are, 
to lowest order in e, 

iDf /f = Vf 



ii9.T 



^Dfff = vf I (i9,_ + l{Af_, + All,)) 
+ {id,, + l{Ali + Afj'^fl 
{id,. + l{At_, + Af_,))f, 



Here, the flux velocity Vf is 



J9 \E\ e 



Vf = 



(5.3) 



(5.4) 



and the covariant time derivative is 



Dl^ = dt + iAll a = 1,2,3,4. 

The same equations hold on 

In the continuum limit, we make the following identifications: 



(5.5) 
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£el _|_ £e2 £e3 i Xe4 

lim = u\ft), lim = v\ft), 

£e3 feA £e2 £el 

lim ^ ^ = u\rt), lim ^ = 

lim Ar_, ^ Al (ft) , lim A- , ^ (rl) , 

/lea 

lim ^ = AUrt), a = 1,2,3,4, 

e^O Vf 



lim J2 = J (5-6) 



on A|,^''. The counterparts to the fields u, v and on sublattice A^^^ are w, z and ^4°, 
respectively. 

The equations of motion, when rewritten in terms of 



u=\ 1, v=\ 1, w=\ 1, z=\ 1, (5.7) 



take the same form as those for four independent massless Dirac spinors in 2+1 space- 
time dimensions, the speed of light being here the fiux velocity. To describe the relativistic 
structure, we need the three gamma matrix 

= {r\\ -rV3, -rV^) = {r\ -\t\ +ir3), 
jp<o) ^ ^M^e(o) ^ (q^ + -.Af)) , (5.8) 

where the r are a new set of Pauli matrices. The long wavelength, low energy limit of Eq. 



( p.llQ when Ji = and for gauge fields close to the fiux phase configuration is 



d^r £, (5.9) 

The dynamics due to fluctuations in the determinant of the W 's has been neglected [Q. 

There are thus four fiavors of Dirac spinors which interact with two independent set of 
SU(2) color gauge fields through the minimal coupling. In the continuum limit, the original 
16 = 2 X (4 + 4) fermionic degrees of freedom per unit cell of A^^^ U A^^^ have become 
16 = 4 X 2 X 2 fields carrying fiavor, color and spinor indices, respectively. On the other 
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hand, the original 16 = 2x4 + 2x4 degrees of freedom residing on the oriented hnks of 
the unit cell have been projected onto the smaller manifold of 4 = 2 + 2 su(2) spatial gauge 
fields. What we have done could have been achieved by expanding the equations of motion 
in reciprocal space around the minima of the mean-field spinon spectrum and keeping only 
the lowest terms in powers of the reciprocal lattice. The gauge fields that are neglected vary 
rapidly compared to the one kept. Choosing a staggered gauge transformation is equivalent 
to choosing the origin in reciprocal space to be the location of one of the degenerate minima. 
The different flavors of the continuum spinous is caused by the degeneracy of the mean-field 
spectrum. 

The local gauge symmetry of Eq. ( |5.9| ) is that of the semi-simple gauge group 
SU(2) xSU(2). Indeed, if we arrange the fermions and gauge fields into the (ji,j2) = 
(1/2, 1/2) reducible representation of SU(2)xSU(2) 



u\ /v\ I Al 



^=1 1, 0=1 1, .1, (5.10) 
then 

£ = Vf ['I' i^ ^' + i^ e] (5.11) 

is invariant under the local gauge transformation 

^(9) ^ U ^(6), -^U A^U^ + (id^ U) U\ 

\ 

W=\ G SU(2) X SU(2). (5.12) 

f/V 

B. Continuum limit in the flux phase when Ji ^ 0. 

The fermions on the even and odd sublattices are coupled by pure fluctuations of the 
degrees of freedom M, Eq. ( |A13|) , on the nearest-neighbor links. These fluctuations are 



of order 0(e) and therefore it is meaningless to separate the contributions from the SU(2) 
color and determinant factors ||5^. Thus, the fluctuating M form a four dimensional real 
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vector space. By repeating the same steps as in the previous subsection, one finds the long 
wavelength, low energy continuum limit of Eq. (|2.11| ) to be given by 



£o + £i + U 



(5.13) 



The first contribution is the one given by Eq. (|5.11|) in the limit Ji = 0. The second 
and third contribution account for the fiuctuations of the M on the nearest-neighbor links. 
There are 16 such degrees of freedom. 

Four of them, ^ . . . , A^|°, transform like scalars under Lorentz transformations. They 
are linearly independent combinations of the original lattice degrees of freedom connecting 
even to odd sites and have been rescaled so has to carry the dimensions of inverse length. 
With the help of the four dimensional representation 

\N{° \ / 
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eof 







eof 







eof 



(5.14) 



the second contribution to Eq. ( p.l3| ) is 



^1 = E tr [AC Ml 

Jl Vf 



(5.15) 



/ / \^ ^ + ^ e + e ^ + e AT^ e 



Notice the mixing of the four Dirac species (or fiavors). Under a local SU(2)xSU(2) 
color gauge transformation U , the A^'s transform like the adjoint representation: 



AC 



U AC U\ n = 1,...,4. 



(5.16) 



As it turns out, the A^'s pick up mean-field expectation values below the fiux phase which 
break the SU(2)xSU(2) symmetry down to the center of SU(2) color. 

The description of the remaining twelve degrees of freedom follows the same line as that 
for the A/''s except that they can be grouped into four families, _B^°, each 
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transforming like a three-vector under Lorentz transformation. With the help of the four 
dimensional representation 



A, 





Aeo 











' 












V. 







r>eo 



(5.17) 



the third contribution to Eq. (|5.13|) is 



= tr 
16 

^ Ji Vf 



A'' A], + B^" Bl + Cl + V^" Vl 



^ / ^ + ^ + e ^ + e 



(5.18) 



The structure of £2 is the same as that for Ci except for the relativistic transformation 
properties of the dynamical fields. Any mean-field expectation value of the A^, . . . , V^'s 
below the flux phase would imply a loss of relativistic covariance. But, as it turns out, 
these fields do not pick up a mean- field expectation value just below the onset of the flux 
phase. We do not expect that they play a relevant role to the issue of the confinement of 
the spinous. 



C. Continuum limit just below {J2/Ji)cv 

Below the onset of the fiux phase, the mean-field solutions of the saddle-point equations 
take the form 

al^O, 4 = 0, 4 = 0, (5.19) 
ReE^O, ImE^O, X 0. 

The six mean- field parameters evolve continuously into their fiux phase values as J2/J1 
approaches from below its critical value. What continuum limit can we expect close to the 
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onset of the flux phase? Strictly speaking, the locations of the minima for the mean-fleld 
excitation spectrum are functions of the parameters 05, E and X, i.e., of J2/ Ji and it would 
appear that any analytical approach is hopeless. 

However, let us expand formally the mean-fleld energy excitations in powers of the lattice 
spacing e by inserting 



ki = koi + ep, / = 1,2,3,4, 



(5.20) 



into the mean-fleld dispersion relation. Here, the lattice spacing e merely plays the role of a 
bookkeeping device. We now assume that 

f lim.^o 7^ \al\ <C, 



3 5, C > 0, 



lim 



lim, 



4 |Re E-Re EA < C, 



1 |Im E-Im Ei\ < C, 



(5.21) 



lim,^o \X\<C, 

so that the formal expansion in powers of e of the mean-fleld excitation energies around any 
of the four points within the Brillouin Zone at which the flux spectrum is gapless, yields the 
relativistic spectra 



(5.22) 



The mean-fleld parameters of Eq. ( [5. 221 ) are taken in the flux phase and 

Ji IXI 



Vf 



J2 iRe E\ e, m = 



2Vf2 



(5.23) 



This simpliflcation suggests that the problem becomes tractable if we can argue that the 
deviations of a\ and E from their flux phase values are negligible compared to the deviations 
of the mean-fleld parameter X from its flux phase value. By inspection of the saddle-point 
solutions, one extracts that on approaching the flux phase from below (see Fig. ^ 
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\al\ oc(l-i^y\ A = 1.52, 



{J2/J: 



1 CI . 



|Re E - Re Efl oc f 1 - ^ ] , /32 = 1.79, 

\ [■hi ■hjcr J 

\lmE-\mE,\^(l- /^a = 1-77, 

The important insight that follows from Eq. (|5.24|) is that jS^ = mm{(3i, ■ ■ ■ , jS/j}. Thus, 
the deviation of X from its flux value is more important than the deviations of a^. Re £^ 
and Im E from their flux values. We will now show how to recover the massive relativistic 
dispersion relation of Eq. ( |5.22| ) from the relativistic theory that was constructed in the 



previous sections when al and J2E are assumed to take their flux phase values whereas JiX 
is non-vanishing. 

In the flux gauge of appendix ^ the sixteen degrees of freedom of Eq. ( |5.14| ) and Eq. 
( p.l7| ) have the mean-fleld values 

4° = 5- = = 5- = 0, /. = 0,1,2, 

Nn = n = l,...,4, (5.25) 

where 

V, = ^ (ia^-W.) = +ie+'5^-, 
v2 

= -^(ia^ + W.) = +ie-'5^-. (5.26) 
v2 

Consequently, at mean-fleld Eq. ( p.l3| ) reduces to 

jC = Co + Ci + £2. (5.27) 

Here, if 
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w 

V 



(5.28) 



then 



+ Vf g mvf M g, 



4 E tr [AC AC! 



0. 



(5.29) 
(5.30) 
(5.31) 



The matrix mM plays the role of a mass matrix for the g's. It is given in a 4x4 dimensional 
representation by 



M 



1 
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(5.32) 



It can be checked that the matrix M, which, in its full glory, is a 16 by 16 matrix, has the 
alternating eigenvalues ±1. Hence, the dispersion relation Eq. ( 5.22| ) is reproduced by Eq. 



(|]23). The global SU(2)xSU(2) color symmetry is broken by Eq. (|5^) all the way down to 



Z2. Again, the generation of the masses for the four Dirac species m, w, v and z is equivalent 
to the breaking of all the continuous color symmetries. 

Finally, by using a linear parametrization of the scalar gauge fields fluctuations, the non- 
linear parametrization of the gauge fluctuations on the next-nearest-neighbor links of Eq. 
( p.l| ) and a linear parametrization of the nearest-neighbor links fluctuations, we find from 
Eq. ( |5.13| ) and Eq. ( |5.3U| ) the intermediate effective theory 



(5.33) 
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D. Symmetries of the intermediate effective theory in the flux phase. 



In the past subsections we showed that it is possible to construct an effective continuum 
theory for the slow modes of the system near the flux regime. The effective theory that we 
constructed has all the right ingredients for the scenario of section |T| to work. In particular, 
below the flux phase, all the continuous gauge symmetries are broken down to a discrete 
Zi and the fermions are all massive. A loop expansion in the gauge coupling constant is 
well defined order by order and yields a local action for the gauge fields. The only low 
energy degrees of freedom left are the adjoint bosonic matter fields (which have acquired 
dynamics through the fermions) and the gauge fields. Hence, the effective action for these 



degrees of freedom will have the form indicated in section |T|. We will present elsewhere the 
calculation of this effective bosonic action . By construction, it will have the two crucial 
features that all continuous symmetries are broken and that the only remaining symmetry 
is discrete namely Z2. In view of the arguments given at the beginning of section p| , the 
spinon quantum numbers are present in the spectrum of the theory provided the effective 
scales for the kinetic energies of the gauge and matter fields are large enough. By varying 
the ratio J2I J\i the parameters /? and k in Eq. ( |3.3| ) of the effective theory change. However, 
the mean-field results alone do not guarantee that it is possible to move from the confining 
regime to the Anderson- Higgs phases of the effective theory. If this can be done, then we have 
a mechanism which does not break parity and time-reversal and allows for the presence of 
the color spinon quantum number in the spectrum of the theory. The only other mechanism 
that we are aware of has its origin in the chiral spin liquid for which the effective continuum 
limit contains a Chern-Simons term fl^. The presence of the Chern-Simons term in the 
bosonized action removes the interactions responsible otherwise for the confinement of the 



fundamental color charge |56 



Let us take the limit J2I J\ 00 for which the Heisenberg model reduces to two inde- 
pendent, nearest-neighbor antiferromagnetic Heisenberg models defined on the even and odd 
sublattices, respectively. Although the exact ground state of the nearest-neighbor, spin-1/2. 
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antiferromagnetic Heisenberg model is not known, there are strong evidences that it carries 



Neel order ^7j. We therefore expect to find degenerate ground states carrying independent 
Neel order on the even and odd sublattices and our mean-field theory must be unstable 
||58| to the Neel ordering. The natural question to ask is whether Eq. ( |5.9| ) can account 
for this instability or if we need to include fiuctuations with higher energy and/or shorter 
wavelength. Since Neel ordering requires the breaking of a global symmetry, let us see what 
symmetries beside the local gauge symmetry Eq. (|5.9|) possesses. 

When J2/J1 00, there exists an additional global U(4) (fiavor) symmetry under the 
transformation 







/u\ 


w 




w 




^ V 




V 




V 









\/V E U(4) = U(l) X SU(4). 



(5.34) 



A possible manifestation of an instability of the fiux phase is the dynamical breaking of the 
global fiavor symmetry as a result of the gauge fiuctuations. We have listed in table | of 
appendix ^ all 15 lattice Hamiltonians Ha expressed in terms of the spinon operators defined 
by Eq. ( p.3| ) which yield all 15 generators qTa q breaking the SU(4) fiavor symmetry. One 
sees that the 9 generators Ti, T3, T4, Tg, T7, Tg, Tn, T13 and T15 are all related to different 
types of valence bond ordering which do not break the spin symmetry . The remaining 6 
generators are related to spin densities and spin currents. For example, q T2 q, measures the 
Neel ordering on the unit cell of the even sites whereas q T4 q measures the Neel ordering 
on the odd sites. Finally, the generator of U(l) is also related to valence bond ordering. We 
thus conclude that Eq. ( ^.9|) has enough symmetries to account for the instabihty towards 
Neel ordering. In fact our effective action could be unstable to many other types of ordering 
including Peierls ordering. 

It is an open problem to show from Eq. (|5.9| ) that the flavor symmetry is broken dynam- 
ically [0 . The difficulty comes from the non-perturbative nature of such a mechanism. For 
example, one can formally integrate the fermions in Eq. (|5.9|) . To one loop, the resulting 
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action for the gauge fields is finite but non-local However, higher loop contributions are 
infrared divergent due to the masslessness of the fermions. The infrared problem can be by- 
passed if one chooses a dimensionless expansion parameter such as 1 /Np with Nf the flavor 
number, instead of the dimensionfuU gauge coupling. Within this approximation, the lowest 
order in l/Np already predicts spontaneous symmetry breaking of SU(4) [^] although this 
result has to be taken with caution When (J2/Ji)cr < -h/ Ji < oo, the effective theory 
Eq. (|5.13|) has only a U(l) global flavor symmetry. For a small number of flavors (such 
as 4 = 2 X 2 in our case), it is likely that the instability of the flux phase with respect to 
the quantum fluctuations at {J2/ Ji)cr is preempted by a flrst order transition from a state 
carrying long range spin-spin correlations to a spin liquid. In any case, Eq. (|5.13|) seems rich 
enough to describe the phase diagram of the Ji — J2 spin-1/2 antiferromagnetic Heisenberg 



model as inferred from a bosonic representation of a Sp(N) generalization of the model p6 



VI. THE ONE DIMENSIONAL CASE. 

In the previous section we have discussed a scenario for separation of spin and charge 
(or the occurence of unusual quantum numbers for excitations) to take place in two dimen- 
sional systems. In this section we want to see what picture this scenario yields about the 
one dimensional Heisenberg chain. Our purpose is twofold. Firstly, the one dimensional 
Heisenberg chain is known to exhibit separation of spin and charge in the sense that it has 
in its excitation spectrum gapless excitations with spin-1/2. These excitations are nowa- 
days referred to as the spinous of the one dimensional chain. Secondly, we want to make 
a comparison with the two dimensional problem. In particular, we want to see what is the 
relation, if any, between the gapless spinous of the one dimensional chain and the spinous 
of the gauge picture. 

The one dimensional Heisenberg chain has been extensively studied by many methods 
and much is known about the physics of this system. The Bethe Ansatz yields the ground 
state and the excitations above it. Conformal fleld theory gives the critical exponents of the 
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correlation functions. The excitation spectrum thus found consists of gapless integer spin 
excitations ("spin waves") and spin-1/2 gapless excitations ("spinous"). 

The integer spin excitations are the expected spin flips. The gapless spin-1/2 excitations 
are the ones we are interested in. The simplest way to see what the "spinous" are is by 
means of the Jordan- Wigner PD| transformation of the spin-1/2 operators: 



c\ = SI e ^J<' J ^ . (6.1) 

From Eq. (|2.17| ), we recognize in the raising and lowering spin-1/2 operators 5*,^ the color 
singlet "baryon" operators 6- and b- , respectively: 

= Sl + iS^ = bj, 
Sr^Sl-iSf = b,, (6.2) 

whereas the "mesonic" operator ml is related to 

S^S. = ^ + = \m^. (6.3) 

The fermion operators q defined by the Jordan- Wigner transformation can be viewed as 
spinless fermions in the sense that they do not carry an explicit spin index. Naturally, the 
expression of Eq. (|6.1| ) depends on the choice of the spin orientation. In this sense the c's do 
transform under the global SU(2) of spin. But, and this is what matters for a comparison 
with the gauge picture, they are singlets under the local SU(2) (color) and consequently 
gauge invariant. Thus, these are allowed excitations which are compatible with confinement 
of the color degrees of freedom. Another important feature of the "spinon" operators of 
Eq. ( |6.1| ) is the fact that they disturb the boundary conditions and, hence, are topological 
excitations. These operators are non-local and gauge invariant. In contrast, the spinous 
of the gauge picture are local but not invariant under local SU(2) (color) transformations. 
In addition, the spinous of the gauge picture, are absent from the physical spectrum if the 
system is realized in a confining state. Clearly, the Jordan- Wigner spinous are not. 
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To sum up, the fermion q, when expressed in terms of the baryons bi and mesons rrii, 
give a non-local and non-linear realization of an operator transforming like the fundamental 
representation of SU(2) under spin rotations and transforming like a color singlet under 
gauge transformations. This is very suggestive of how confinement precludes any simple 
description of the spin-1/2 low energy sector of the one dimensional Heisenberg model in 
terms of color carrying spinous. 

Let us see what the SU(2) gauge approach tells us about the one dimensional case. In 
particular we would like to see how the known physics of this system is recovered in the 
gauge field approach. It is straightforward to apply the methods of the previous sections 
to the one dimensional case. In a separate publication we do so |^. Here, we will discuss 
several qualitative differences that arise when this picture is applied to the one dimensional 
case. 

Along the lines of section [IV C|, we try the mean-field Ansatz 



Ao^t = A„ (6.4) 

f -Xa^ if j = i + l, 

—Re E + Im E a"^ if j = i + 2, i even, 

-Re E -ImE ifj = i + 2,i odd, 

otherwise. 
It describes a spin liquid if and only if 

Ao = ^4 a\ (6.5) 

This Ansatz satisfies all of the requirements: a) it is translationally invariant (up to gauge 
transformations), b) it respects all point group symmetries and time reversal invariance (up 
to gauge transformations) and c) it achieves a complete breaking of the local continuous 
gauge symmetry through the operators 

P: = W^.(,+i)W^(.+i)(.+2)W^W2)' (6-6) 

Pi+l = W^(^+l)(^+2)W^(^+2)(^+3)W^(^+l)(^+3)' ^VeU. 
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In particular, just as in the two dimensional case, it leaves a discrete symmetry unbroken. 

However, discrete gauge symmetries play a very different role in one and two space 
dimensions. This is so since the lower (space-time) critical dimension for gauge theories 
with a discrete symmetry group (such as in the case of interest for our problem) is three. 
Hence, in the case of the one dimensional quantum antiferromagnet, which lives in two space- 
time dimensions, the effective Z2 gauge theory along the topmost line of the phase diagram 
in Fig. |l| is below the lower critical dimension and thus is always in a confining phase. No 
phase transition is encountered along the line (/5, k. A) = (/?, 00, 00) and there can only be 
one phase in the phase diagram, the confinement-Higgs phase. In this phase, the spinon 
carrying the fundamental color charge is confined and is not part of the spectrum of finite 
energy states. Clearly, the spin-1/2 excitations of the one dimensional quantum Heisenberg 
antiferromagnet are not the spinous of the gauge picture, as has been claimed rather loosely 
in the literature of the subject. As we stressed above, the fermion states described by the 
Jordan- Wigner operators are manifestly gauge invariant under the local SU(2) symmetry 
and are unaffected by confinement. 

VII. CONCLUSION. 

We have studied the issue of spin and charge separation for strongly correlated electronic 
systems in two spatial dimensions by studying the frustrated Heisenberg model for spin- 
1/2 in a representation of the spin degrees of freedom in terms of fermions (spinous) and 
gauge fields. Our treatment preserves explicitly the existing SU(2) local gauge symmetry 
which allows us to treat on the same footing the Affleck-Marston and the Anderson order 
parameters. We have given a set of sufficient conditions for the system to have excitations 
in which the spinous carry the separated spin quantum numbers, i.e., spin-1/2 degrees of 
freedom, inspite of the existence of strong gauge fluctuations. 

We solved the saddle-point equations for an isotropic mean-field Ansatz first proposed 
by Wen and depending on four real parameters. There exists a second order phase transition 
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to a new flux phase which is triggered by the competition between the nearest Ji and next- 
nearest-neighbor J2 antiferromagnetic couphngs of the model. The flux phase is chosen for 
strong enough J2. Below the flux phase, a gap for the spinous opens and reaches a maxima 
for J2/ Ji ~ 0.7. Although no level crossing occurs for the mean-fleld ground states below the 
flux phase, there is a level crossing for the excited states when J2/ Ji ^ 0.7. On approaching 
the pure nearest-neighbor limit, the features of a tight-binding Fermi surface at half-flUing 
emerge. The phase < J2/J1 < {J2/ Ji)cr, the s-RVB phase, is characterized by an energy 
gap and the breaking of all continuous gauge symmetries. 

We derived a continuum theory for the soft degrees of freedom in the vicinity of the flux 
phase which describes bosonic gauge and matter flelds interacting with the spinous. In the 
flux phase, the continuum theory is likely to be unstable towards formation of long range 
antiferromagnetic order or dimer ordering. A more detailed calculation is needed to decide 
which ordering is chosen by the system but, in any case, deconflnement of the spinous is not 
expected to take place. We showed that below the flux phase, the only important remaining 
soft modes are those of a Z2 gauge theory in three space-time dimensions due to the break- 
ing of all continuous gauge symmetries by the s-RVB Ansatz. Deconflnement of the gauge 
spinous is then allowed in the phase diagram of a generic theory with the same symmetry 
attributes but in which the gauge and matter fields have independent kinetic energy scales. 
In our case, the calculation of the ratio of the effective energy scales for the gauge and matter 
flelds through the integration of the spinous is needed to decide whether the system is in 
the confining or deconfining regime. Consequently, although deconflnement of the spinons 
is allowed it is not guaranteed since it depends on the values of microscopic parameters. 
This result stands in sharp contrast to the analysis of the one dimensional quantum anti- 
ferromagnet which predicts that the gauge spinons are always conflned, irrespectively of the 
value of the effective energy scales for the gauge and matter flelds, although the system still 
exhibits separation of spin and charge. The mechanism for separation of spin and charge in 
one dimension is topological and it is unrelated to the issue of confinement. 
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APPENDIX A: THE FLUX PHASE 

We have found in section ^ that above a critical value of the coupling constant J2/ Ji, 
the s-RVB phase collapses to the case 

as = 0, Re E = \m. E, X = 0, (Al) 

i.e., the spectrum takes the simple form 



l^fcl = J2 |Re -E"! \J cos'^ kx cos'^ ky + sin /c^sin ky. (A2) 

The spectrum is thus gapless and linearization around the four minima of Eq. ([4.14|) yields 
a massless relativistic dispersion relation. 

Another characterization of this phase can be made. We introduce the orthogonal and 
normalized basis 

= >n = -^. W^ = -.'. (A3) 

of the Lie Algebra su(2). This basis satisfies the su(2) Algebra with 

W-W+ = iWs, W+W3 = iW-, WsW- = iW+, (A4) 
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and can be used to rewrite 

' along the links x and y, 

W^j^ = I \E\W+ along the link x+, (A5) 

^ along the link x_, 

whenever J2 is larger than its critical value. This regime is characterized uniquely by the 
most elementary path ordered operator along a closed path: 

= W_W^Wlwl = -cr°. (A6) 

Here, is the neutral element of SU(2). Eq. ( [A^ ) tells us that there is a mean-field fiux of 



71 per elementary plaquette, a situation commonly referred to as the fiux phase [|14| , p!5 . 

As of now, the mean-field Ansatz had an explicit translational invariance. In order to 
describe small fiuctuations of the link and time-like su(2) degrees of freedom around their 
mean-field value in the fiux phase, it turns out to be convenient to choose a gauge in which 
the color structure on the next-nearest-neighbor links becomes trivial. The price to be payed 
is a reduction of the translational symmetry. First we need some notation. 

We introduce eight sublattices (see Fig. H) by partitioning first the square lattice A, 
whose sites i E are labelled by pairs of integer, into two interpenetrating sublattices Ag 
and A^: 

Ag = {ieA\ii + i2 = mod 2 } , 

A^ = {ie A\ii+i2 = l mod 2 } . (AT) 

The even and odd sublattices are then partitioned into four sublattices A.'f'\ A'^\ cl = 
1, 2, 3, 4 by choosing 

A« = {ieA,\3jeZ',t = n2x.+j22x+}, 

A« = {2GAJ3jgA«,z = j + x}, (A8) 

and using translations to define 
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A(^) = {ieA,|3ieA«,i = j + x_}, 

A(^) = {ieA,|3jeAf),i = j-x_}, (A9) 

and similarly for the odd sublattices and A^'^^. From now on, i will always be taken 

to belong to A^^). 

The fermionic doublets will be renamed 

"^i+x- ~^ fi ' '4'i+x+x- ~^ fi -I 

'ri+a;_+x+ J I 1 yi-\-x-\-X--\-x+ Ji ' 

A+x^ - ft, A+x+x^ - ft- (AlO) 

In the same way, there will be 2 x 4 = 8 independent time-hke gauge fields: 

~^ -^ti-i " ' 1 '^6{i+x+x+) ~^ > (AH) 

2 X 8 = 16 next- nearest- neighbor link degrees of freedom: 

- • • • (A12) 

■ ■ ■ ' ^{i+x+x+){i+x+2x+) ~^ Qx+i, 

and 2 X 8 = 16 nearest-neighbor hnk degrees of freedom: 

W,^+,)^Mll,... (A13) 

■ ■ ■ W — > 

' ^'^ {i+x+x+){i+x+x++y) ^^^yi ■ 

We are now ready to perform the staggered gauge transformation 
ft ^ sgn(z) (+<7°) ft\ ft ^ sgn(z) (+<7°) 
/f ^ sgn(^) (-iW-_)/f , /f ^ sgn(z) {-iW.)ft. 
ft ^ sgn(i) {-m)ft, ft ^ sgn(i) (-11^3)/^, 
ff ^ sgn(0 (+iiy+)/f , - sgn(z) (+iiy+)/f , 



(A14) 
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for all i G A^^\ Here, sgn(i) is the sign function which distinguishes between even and 
odd sites of sublattice A^^^. Under this gauge transformation, the mean- field Ansatz for the 
nearest-neighbor links takes the simple form 



Qx-i 


= +i\E\ 


<y\ 




= 


a° 


Qx-i 


= +i\E\ 


a', 


Qx+i 


= 


a° 


Qx-i 


= +i\E\ 


a\ 




= 


a° 


Qx-i 


= +i\E\ 




Qx+i 


= 





and similarly on One easily verifies that Eq. ( [A6D still holds. 



APPENDIX B: CONTINUUM LIMIT FOR SOME LATTICE HAMILTONIANS. 

We have collected in this appendix all the 15 lattice Hamiltonians which yield the 15 
generators of the global SU(4) symmetry of the Lagrangian density Eq. ( ^.9|) . They are 
displayed in table | below. 
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FIGURES 

FIG. 1. Phase diagram for a gauge-Higgs lattice theory in 2+1 dimensional space-time. /3 is 

the energy scale for the gauge fields, k that for the matter (Higgs). The amplitude of the Higgs is 
frozen (A = oo). (a) shows the case of the matter (Higgs) in the fundamental representation, (b) 
shows the case of the matter (Higgs) in the adjoint representation. 

FIG. 2. The solutions ai, Re Im and X as a function of the dimensionless ratio J2/J1 
of the saddle-point equations for Wen's Ansatz. 

FIG. 3. The upper branch of the MF single-particle excitations (spinous) in the first 

quadrant of the Brillouin Zone for: (a) J2/J1 = 1.4, (b) J2/J1 = 0.9, (c) J2/J1 = 0.5, (d) J2/J1 = 0. 
The minima move from the edges (a) and (b) to the diagonal (c) and (d) . 

FIG. 4. Phase diagram of Wen's mean-field Ansatz. There exists a second order phase 
transition into a phase, the flux phase, without any gap to the mean-field excitations. Below 
the threshold (J2/Ji)cr ~ 1.342, the mean-field excitations have a gap, except at the BZA point 
J2/J1 = 0. The critical point {J2/Ji)cr is infrared unstable as indicated by the RG flow of J2/Ji- 
The location of the minima in the s-RVB phase, < J2/J1 < {J2/ Ji)ct^ jumps discontinuously for 
J2/J1 Ki 0.7, signaling a level crossing of the excitations. The values of the gap are plotted in units 
of Ji. 

FIG. 5. The scaling exponents /32, /^s and /34 of the s-RVB parameters a\, Re E, Im E 
and X as a function of J2I J\ can be read from a log-log plot. We have used {J2/ Ji)cr = 1.342 and 
Re Ei = lmEi = 0.339. 

FIG. 6. Partitioning of the square lattice A into eight sublattices. A sublattice is labelled by 
the letter e or o for even or odd respectively and by an integer between 1 and 4. 



54 



TABLES 

TABLE I. The 15 bilinear s f (fx q Ta q are given together with their lattice counterparts. 
Here, Ta, a = 1, • • • , 15 are the generators of SU(4) flavor. The spinons s were introduced in section 
II and are relabelled according to the lattice partitioning of appendix A. The u, w, v and z are 
defined in section V. 

a Lattice Hamiltonian Ha J <fix g g 

1 EisA« [<'^^°«f + Si*^^°«f - H-C] + / d^x [uv + vu] 

2 +EieAe -^f^'^^^f +H.C.] -if d^x[uv-vu] 

3 -iEieAe +sfa^sf -H.C.] +/ d^x[uu-vv] 

4 -iE^sAj [sf'-'sf +sf^a°sf -H.C.] +/ d^x[v 

^zGA^ L * I — 2a;_|_ — 2a:_ i % — 2xj^ ' « i— 2x_|_ z z — 2x_|_ 

-^ft^o^oi -sf^o^f +sftaO<i +sfVo<2^^_ -H.C.] +/ <fx[uw + wu] 



10 



13 



Z_-/zeA| L * * * I ' I %—2x— I 



\wz + ^u; 



-sl^t^S^ol _^e3t^3^o2 +^^41^3^01 + ^^-^f ^3^o2 ^ + H.C.] -i / (fix[ 



uw — wu\ 



- J d'^x [uz + zu] 



-5^=3+^3503 _5e3t 3^04 +s«4t^35C,3 + 5^4! 3^04 + H.C.] -if d'^xlv 

I I I X ' I l — 2x— ' 2 I ' J J I- 



+sfa^sf - sfaOsf + ^'o^st,,_ - sfa<^sf - H.C.] + / <Px [ 



wv + vwl 



/ j-iez L i I ' I i — lx_ % % % % 

+S'r^"''cr3s°3 _s'=3t^35°4 +s'r4t^35o3 _ s'r4t^35o4 + H.C.I -if d^a; te^, _ ^^1 
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+ j (fix [vz + zv\ 



14 



15 



■E 



. o- o- +s'^^^a'^s°'^„. +H.C. 

2a:_i_— 2a;_ ' i z— 2a;_i_ ' 



-^^■=Ae [s"'V05=3 +5'flt^03e4 _5e2t O^eS _ ^=21 0^e4 _ ^ _ H.C." 
2 Z_/2(=A^ z z I I z z z^ 



— i J (fix \pz — zv\ 
+ J" d^o: [i;?; — zz\ 
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